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Abstract 

The maximum number of non-crossing straight-line perfect matchings that a set of n 
points in the plane can have is known to be 0(10.0438") and e*(3”). The lower bound, 
due to Garcia, Noy, and Tejel (2000) is attained by the double chain, which has 0 ( 3 "n‘^O)) 
such matchings. We reprove this bound in a simplified way that uses the novel notion of 
down-free matching, and apply this approach on several other constructions. As a result, we 
improve the lower bound. First we show that double zigzag chain with n points has 0*(A") 
such matchings with A « 3.0532. Next we analyze further generalizations of double zigzag 
chains - double r-chains. The best choice of parameters leads to a construction with 0*(j^") 
matchings, with v « 3.0930. The derivation of this bound requires an analysis of a coupled 
dynamic-programming recursion between two infinite vectors. 
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1 Introduction 

Background. A non-crossing straight-line matching of a finite planar point set is a graph 
whose vertices are the given points, whose edges are realized by pairwise non-crossing straight 
segments, and where every vertex has degree at most 1. In what follows, such matchings will be 
simply called matchings. A matching is perfect if every point is matched - that is, has degree 1. 
Throughout the paper, all point sets are assumed to be in general position in the sense that no 
three points lie on a line. 

In this paper we deal with bounds on the number of perfect matchings that a set of size n can 
have. This question arises in a broader context. Non-crossing straight-line matchings, either 
perfect or not necessarily perfect, are just two kinds of geometric plane graphs, others being 
triangulations, spanning trees, connected graphs, etc. A web page of Adam SheffeiQ maintains 
the best up-to-date bounds on the maximum number of geometric plane graphs of several kinds. 

First we recall that for the minimum number of perfect matchings that n points in general 
position can have, the exact solution is known. Garcia, Noy, and Tejel [6] proved the number of 
perfect matching is minimized on point sets in convex position. It is well-known that the number 
of perfect matchings is then 0^/2, where Ck = = 0(4^//c^/^) is the A:-th Catalan number. 

The minimum number 0^/2 of perfect matchings is in fact attained only for point sets in convex 
position, with the exception of one configuration of six points [2] . 

Regarding the maximum number of perfect matchings that a point set of size n can have, 
only asymptotic bounds are known. The best upper bound to date, 0(10.0438”), was proved 
by Sharir and Welzl [9]. The best previous lower bound was given by Garcia, Noy, and Tejel in 
the above-mentioned paper |^. They showed that for the so-called double chain with n points 
(to be denoted by DC^,) the following holds: 

Theorem 1 ([6l Theorem 2.1]). The number of perfect matchings of the double chain with n 
points is 0 (3”n‘^(^)). 

Actually, it follows from their proof that this number is n(3”/n^) and 0(3”/n^). In Sec¬ 
tion we shall sketch this proof, and also determine the polynomial factor more precisely. 

The double chain was used in [6] not only for improving the lower bounds on the maximum 
number of perfect matchings, but also for some other kinds of geometric graphs: triangulations, 
spanning trees and polygonizations. It was believed by some researchers in the field that it 
might give the true upper bound at least for some of these kinds. However, in 2006, Aichholzer, 
Hackl, Huemer, Hurtado, Krasser, and Vogtenhuber [T] introduced a new construction, the 
double zigzag chain with n points (to be denoted by DZZC„), see Figure]^ below. They proved 
that DZZC„ improves the lower bound for the number of triangulations: it is 0*(8”) for DC^ 
and 0*(8.48”) for DZZC„,. (The notations O*, Q,*, and 0* correspond to the usual 0-, 12- and 
0-notations, but with polynomial factors omitted.) To our knowledge, the number of geometric 
graphs of other kinds mentioned above for DZZCn was not found. 

^Numbers of Plane Graphs. http://www.cs.tau.ac.il/~sheffera/counting/PlaneGraphs.html 
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In this paper we determine asymptotically the number of perfect matchings for DZZC^ and 
its further generalizations, improving the existing lower bound. 

Our results. In Section we will first show that DZZC,i has asymptotically more perfect 
matchings than DCn: 

Theorem 2. The number of perfect matchings of the double zigzag chain with n points is 0*(A”’), 
where A = ^J{^/m + 9)/2 ^ 3.0532. 

In Sections andwe will present a generalization of DCn, which comes in two variations: 
double r-chains without corners and double r-chains with corners, see Figures and below. 
Our best results for these constructions are as follows: 

Theorem 3. The number of perfect matchings of the double 11-chain without corners with n 
points is where v = V240054 « 3.0840. 

Theorem 4. The number of perfect matchings of the double 8-chain with comers with n points 
is D((zz—e)”-), and where v = ^/(8389~-|-~3\/777T7i7y^ ss 3.0930 and e > 0 is arbitrarily 

small. 

We shall present proofs for all three theorems because they use different techniques. First, 
in Section we introduce the notion of down-free matchings and show in Theorem how one 
can generally reduce the problem of asymptotic enumeration of perfect matchings of a “double 
structure” to that of down-free matchings of the corresponding “single structure”. In the proof 
of Theorem]^ (Section]^, we find a recursion for the number of down-free matchings of the 
zigzag chain, and translate it into a functional equation satisfied by the generating function. We 
solve this equation explicitly, which allows us to find the asymptotic growth rate by looking at 
the smallest singularity of the function. In the proof of Theorem]^ (Section]^ we use matchings 
which possibly have runners - edges with only one endpoint assigned. We define a sequence 
of infinite vectors whose entries are the numbers of down-free matchings of the r-chain of a 
certain size, sorted by the number of runners. These vectors can be computed recursively. We 
reformulate this recursion in term of lattice paths and obtain the desired growth rate with the 
help of a result of Banderier and Flajolet [3]. The proof of Theorem]^ (Section]^ starts similarly, 
but due to technical obstacles, we need two sequences of infinite vectors, defined by a coupled 
recursion. We find that the desired growth constant is determined by the dominant eigenvalue 
of certain 2x2 matrix. 

Notation. We use the following notation and convention. A construction A is a family 
{Xn}n£i for some infinite I C N, where, for fixed n, is a class of point sets of size n with 
certain common properties, for example, a certain order type (or, in some cases: one of several 
order types) and certain restrictions concerning position in the plane with respect to coordinate 
axes. The double chain (DC) mentioned above is one of such constructions. Occasionally we 
will abuse notation and denote by A„ not only such a class, but also any of its representatives. 
If we know that all members of A„ have, for example, the same number of matchings, we can 
speak unambiguously about “the number of matchings of , and so on. 

In what follows, pm(A„) denotes the number of perfect matchings of am(A„), the 
number of all (non-crossing straight-line, but not necessarily perfect) matchings of dfnn(A„), 
the number of down-free matchings of A„. For some constructions it can happen that not all 
representatives of Xn have the same number of (for example) perfect matchings and, thus, 
pm(A„) is not well-defined, but the common asymptotic bound still can be given, which enables 
us to write expressions like pm(A,i) = Q*{gA) in such cases as well. 

For two distinct points p and q, the straight line through p and q will be denoted by £{p, q). 
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A set of points is in downward position (respectively, in upward position) if the points lie 
on the graph of a convex (respectively, concave) function. In particular, three points with 
different x-coordinates are in downward position (respectively, in upward position) if they form 
a counterclockwise (respectively, clockwise) oriented triangle when sorted by x-coordinate. 

A point of X not matched by a matching will be called a free point. 


2 Double chains and double zigzag chains 

In this section we recall the definitions of a double chain and a double zigzag chain, and recall 
how the bound pm(DCn) = &*from Theoremfor was obtained in [6]. 

2.1 One set high above another and “double constructions” 

We start with a notion of a point set being one “high above” another point set. 

Definition. Let P and Q be two point sets in the plane, each of them having distinct x- 
coordinates. We say that P is high above Q if any point of P lies above any line through two 
points of Q and any point of Q lies below any line through two points of P. 

Observation 5. For any two point sets P and Q with distinct x-eoordinates it is possible to put 
a translate of P high above a translate of Q. 

Proof. Assume that Q lies in the interior of a “topless rectangle” R = [a,b] x [c,+oo). Any 
line that passes through two points of Q crosses the lines x = a and x = b in certain points. 
Therefore there exists a number d such that all points of R above y = d are also above any line 
containing two points of Q. A similar argument (with “below” instead “above”) applies for P. 
Thus, if we put translates of P and Q between vertical lines x = a and x = h (where a and b are 
chosen to ht bot sets) and, if needed, translate P upwards, we get eventually a translate of P 
high above a translate of Q. □ 

Let Xn be a construction. A double Xn (denoted by P)X 2 n) is the family of sets obtained 
by placing P, a representative of high above Q, a representative of Xn reflected across a 
horizontal line. An edge between a point of P and a point of Q will be called a PQ-edge. 

2.2 Double chains 

A (single) downward ehain (respectively, upward chain) of size n is a set of n points in downward 
(respectively, downward) position. A downward chain of size n will be denoted by SC^. 

Let n be an even number. A double chain of size n consists of a downward chain of size n/2, 
P = {pi,P 2 , • • • jPn/ 2 }) placed high above an upward chain of size n/2, Q = {qi,q 2 , ■ • • ,(?n/ 2 }- 
See Figurefor an example. A double chain of size n will be denoted by DC^. We assume that 
both sets P and Q are sorted by x-coordinate. 

2.3 Perfect matchings in the double chain 

Theorem was proved in [6] as follows. 

Denote by pm(DC„j) the number of perfect matchings of DC,i that have exactly j PQ-edges 
between the upper and the lower chain. If n/2 — j is odd, then no perfect matching exists, so 
we assume that n/2 — j is even. One can construct a perfect matching with j PQ-edges in 
the following way. First choose any j points of P and j points of Q and connect them by j 
non-intersecting PQ-edges. It is easy to see that there is a unique way to connect the chosen 
points (see also Proposition below). Then, choose any perfect matching of the free points in 
each chain. Alternatively, one can first choose n/2 — j points of P and n/2 — j points of Q, then 
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P 




Figure 1: A double chain of size 22. 


take any matching of P and any matching of Q that uses the chosen points; after that, the free 
points can be matched by PQ-edges in a unique way. Since Q has the same order type as P, it 
follows that ^ 

pm(DC„,,) = (am(SC„/2,,))' = • C^n/2-j)/2'^ , (1) 

where am(SC„/ 2 ,j) denotes the number of matchings of P (or, equivalently, of any set of n/2 
points in convex position) with exactly j free points. Finally, the total number of perfect 
matchings of DC^ is 

pm(DC„)= ■C^n/ 2 -j)/ 2 ) • (2) 

0<j<n/2 \ \ / / 

j=nl2 (mod 2) 

An analysis shows that the dominant term in this sum is the term corresponding to j ~ n/6, 
that is • C'n/e) (n/6 should be rounded in one way or the other). Using the estimates 

ibk) ~ and Ck = 0(4^/fc^/^), which follow from Stirling’s formula, one 

obtains pm(DC„^„/6) = ©(S^/n^), and, therefore, pm(DCn) = n(3”'/n^) and 0(3”/n^). With 
the help of Stirling’s formula, and replacing the sum Q by an integral, one can obtain the more 
precise estimate pm(DC,i) = 3”/n 7/2 . 182/7r3/2 . (1 +o(l)) (we omit the details). 

2.4 Double zigzag chains 

In this section we recall the dehnitions of a (single) zigzag chain SZZC and a double zigzag chain 
DZZC. 

Let P = {pi,P 2 , ■ ■ ■ ,Pn} be a downward chain (SC„) sorted by x-coordinate. For each even i, 
1 < i < n, we move the point pi vertically up, very slightly above the segment pi-ipi^i, so that 
all triples pi-ipipi+i with even i {1 < i < n) are now in upward position, and all other triples of 
points remain in downward position. After this modihcation, the points pi,p 2 , ■ ■ ■ ,Pn are still 
sorted by x-coordinate. A set obtained in this way will be called an even (single) downward zigzag 
chain of size n and denoted by eSZZC,i. If instead of even is we perform this transformation 
for each odd i, 1 < i < n, we obtain an odd (single) downward zigzag chain (oSZZCn). If n is 
even, then eSZZC^ and oSZZC^ are reflections of each other with respect to a vertical line; but 
if n is odd, then eSZZC^ and oSZZCn have different order types, and - as one can verify on 
some small examples - different numbers of (perfect or not necessarily perfect) matchings. See 
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Figure [^for an example. A zigzag chain of size n, denoted by SZZC„, is either an eSZZC„ or an 
oSZZCfi- For both types of SZZC^, we shall derive the same asymptotic bound on the number 
of perfect matchings. 

An even DZZC for odd n An even DZZC for even n 



An odd DZZC for odd n 


An odd DZZC for even n 



Figure 2: A (single) zigzag chain - several cases. 


An upward zigzag chain (of either kind) is a downward zigzag chain reflected across a hor¬ 
izontal line. The construction of a double zigzag chain from zigzag chains is analogous to the 
construction of the double chain from two single chains: A double zigzag chain of (even) size n 
(DZZCn) consists of a downward zigzag chain P = {pi,P2, • • • ,Pn/ 2 } high above an upward 
zigzag chain Q = {gi; 92, • ■ •, Qn/ 2 }- We can combine even and odd zigzag chains in various 
ways, but as mentioned above, this will make no difference for the asymptotic number of perfect 
matchings. See Figure for an example of double zigzag chain obtained from two even zigzag 
chains of odd size. 


P 





Figure 3: A double zigzag chain of size n = 22. 


3 Down-free matchings and perfect matchings 

3.1 Down-free matchings 

Suppose that we want to adapt the argument that was used for estimating pm(DCn) for the case 
of pm(DZZCn) (of any kind). That is: for hxed j (such that n/2 — j is even) we want to choose j 
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PQ-edges, and to complete this matching to a perfect matching by choosing edges that connect 
free points of the same chain in all possible ways. One can hope for improvement since the 
number of perfect matchings in SZZC^ (of any kind) is with n = \/2 + 2\/2 ss 2.1974, 

in contrast to 0*(2”') for SC„. (This bound for SZZCn was proven in [T] for a slightly different 
construction, the so called double circle. The order type of SZZC is different from that of a 
double circle only in one triple of points; it is easy to show that they have the same asymptotic 
number of perfect matchings.) However, in comparison with the case of DC^, here we have 
less freedom and no uniformity in constructing the matchings inside P and Q, once PQ-edges 
are chosen. Indeed, the j chosen PQ-edges may block visibility between certain pairs of free 
points from P or from Q. Moreover, for different choices of j PQ-edges, we have in general 
different numbers of ways to complete them to a perfect matching of DZZC„. This follows from 
the fact that sets of points that remain free after choosing j PQ-edges have in general various 
order types, and, so, it seems hopeless to enumerate them in this way. On the other hand, if we 
first choose (n — 2j)/4 edges between two points of P and (n — 2j)/4 edges between two points 
of Q, then - as we prove below in Proposition!^ - there is at most one way to complete such a 
matching to a perfect matching of DZZCn. More precisely, if the free points of P “see” all free 
points of Q, there is exactly one way of complete a matching to a perfect one, otherwise it is 
impossible. Next we define a property of matchings which - for two sets being one high above 
another - ensures the desired visibility of free points. 

Definition. Let P be a set of points with distinct x-coordinates. Let M be a matching of P. 
AI is a down-free matching if for each unmatched point p G P, no edge of M has a point directly 
below p. In other words: for each free point p ^ P, the vertical ray going down from p, does not 
cross any edge of M. Similarly, one defines an up-free matching. 

Proposition 6. Let P and Q be two point sets in general position with distinct x-coordinates 
such that P is high above Q. Let Mp be a matching of P and Mq a matching of Q, such that 
Mp and Mq have the same number of free points. 

1. If Mp is a down-free matching and Mq is an up-free matching, then Mp L) Mq can be 
completed to a perfect matching of P U Q in a unique way. 

2. If Mp is not down-free or Mq is not up-free, then it is impossible to complete MpU Mq 
to a perfect matching of PL) Q. 

Proof. We assume again that the points P = {pi,P 2 ,P 3 , • ■ ■ } and Q = {qi, q 2 , Qs, ■ ■ ■} are sorted 
by x-coordinate. 

First we observe that for any Pa,P/3 £ P, q-y^QS £ Q, the points Pa,P^,qy,qs are in convex 
position. Indeed, if, for example, qs G comr{pa,pp,qy), then the points pa and pp would lie on 
different sides of the line qs), and this contradicts P being high above Q. 

Now we show that for a < /? and y < 5 the points Pa,Pi3,qs,qy lie on the boundary of 
their convex hull in this clockwise order, see Figure |^a) for an illustration. Since P lies high 
above Q, the points of Q lie below I{pa,Pp) and thus the points pa and pp lie on the convex hull 
consecutively and in this clockwise order. Similarly qs and lie on the convex hull consecutively 
and in this clockwise order. This implies the claim. 

I. Assume that Mp is down-free and Mq is up-free. 

Let Pax , pa 2 , ■ ■ ■,Pam b® t^e free points of P and let q.yxNy 2 ^ ■ ■ ■ ^ q^m be the free points 
of Q, sorted from left to right. We complete Mp L) Mq to a perfect matching of P U Q 
by connecting with q-/. for i = 1,2,... ,m. By the just-proven claim about the cyclic 
order of Pa,Pp, qs, qy, these new PQ-edges do not cross each other. Moreover, they do not 
cross the edges of Mp and of Mq. Indeed, assume that an edge Paqy = PaNyx crosses an 
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edge e € Mp. Consider the angular sector F bounded by the downward vertical ray ri 
with the origin pa and the ray r 2 from pa through see Figure [^b). The edge e crosses 
the ray r 2 by assumption and does not cross the ray ri, because the matching Mp is 
down-free. Therefore, one of the endpoints of e, say px, lies in the interior of F. However, 
this is impossible because in such a case q-y lies above the line i{pa,Px)- 

Finally, we need to show that this is the only way to complete Mp U Mq to a perfect 
matching of P L) Q. Indeed, for any other possibility to match the free points we would 
have a pair of edges PaQs and p^q-y with a < fi, y < 5. However, it follows from the claim 
about the cyclic order of Pa,Pi 3 ,q 5 , q-f that such edges necessarily cross. 

2. Assume without loss of generality that Mp is not down-free. Then there is a free point p^ in 
Alp so that the vertical downward ray from pjs crosses an edge PaP-f, with a < j3 < y. See 
Figurej^c) for an illustration. The set Q must lie below i{pa-,pp)^ ^{Pa,P'y), and i{pp,p-y)- 
There is no way to connect pjs to a point q G Q without crossing the edge paP-y ■ D 

3.2 Down-free matchings of X and perfect matchings of double X 

In the following theorem we show how asymptotic bounds on dfm for a structure X imply those 
on pm for the corresponding double structure PX. 

Theorem 7. Let X be a eonstruction so that dfm(X„) = 0*(A’^). Then for the double structure 
DA we have pm(DA„) = 0*(A’^). 

More precisely: //dfm(A„) = 0(A”/n“), then pm(DA„) = D(A"'/n^“'''^) and 0(A’^/n^"). 

Proof. Denote by dfmj(A„/ 2 ) the number of down-free matchings of A „/2 with exactly j free 
points, for 0 < j < n/2, and let pj = dfmj(A„/ 2 )/dfm(A„/ 2 )- Then we have Ylo<j<n/ 2 Pj ~ 
which implies < J 2 o<j<nl 2 P'j < 1- Now 

pm(DA„) = pmj{PXn)= Y = dfm(A„/ 2 )^ • Y Pp 

0<j<n/2 0<j<n/2 0<j<n/2 

which implies the claim immediately. □ 

As the first application of Theorem we show how one can reprove Theorem [^without need 
to determine the dominant term in Equation Q. We use the following well-known fact. 

Proposition 8 ([8] A001006). The number of all matchings in a set of n points in convex 
position is the nth Motzkin number Mn. Asymptotically, M^ = 0(3"'/re^/^). 

Moreover, every matching of a downward chain is obviously down-free. Therefore, Theoremj^ 
with A = 3 and a = 3/2 gives immediately pm(DC„) = D(3"'/n^) and 0(3"/n^). 

In the next sections we use Theorem [7] for estimating pm for other constructions. 
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4 Zigzag chains 

By Theorem asymptotic bounds on dfm(SZZC„) imply those on pm(DZZCn). Thus, we 
analyze the number of down-free matchings of SZZC„. We defined above two kinds of double 
chains: even and odd. We introduce three generating functions depending on the kind of chain 
and on the parity of n: 

1. A{x) = where = dfm(eSZZC 2 fc-i-i); 

2. B{x) = where bk = dfm(oSZZC 2 fc+i); 

3. C{x) = Ylk>o^kX^, where = dfm(eSZZC 2 A:) = dfm(oSZZC 2 fc)- 

We find recursive relationships between the coefficients of these functions. 

Recursion for a^. For every k > 0 we have the following cases (Figure]^. 

1. Pi is not matched. This contributes Ck matchings. 

2. Pi is matched to P 2 i-i-i with 2 < i < k. This contributes Yl, 2 <i<k h-iCk-i matchings. 

3. Pi is matched to p 2 i with 1 < i < k, p 2 i-i and P 2 i-i-i are not matched to each other. This 

contributes Yli<i<k matchings. 

4. Pi is matched to p 2 i with 2 < i < k, P 2 i-i and P 2 i-i-i are matched to each other. This 
contributes ^ 2 <i<k ^i- 2 Ck-i matchings. 

5. Pi is matched to ps. Then p 2 must be matched to some point P 2 i+i with 2 < i < k. This 
contributes Yli 2 <i<k ^i-^Ck-i matchings. 

6. Pi is matched to ps, p 2 is matched to p 2 i with 2 < i < k, and P 2 i-i and P 2 i+i are not 
matched to each other. This contributes Yli 2 <i<k Ci- 2 <ik-i matchings. 

7. Pi is matched to pa, p 2 is matched to p 2 i with 3 < i < k, and P 2 i-i and P 2 i-i-i are matched 

to each other. This contributes Yl 3 <i<k matchings. 

Thus we obtain 

Oik — Ck~\~ ^ ^ bi—\Ck—i~\~ ^ ^ Cj_ifl/c_jT2 ^ ^ hi—2Ck—i~\~ ^ ^ Cj—20-^—iT ^ ^ bi—^Ck—i- (3) 

2<i<k l<i<k 2<i<k 2<i<k 3<i<k 

Recursion for bk- For every A: > 0 we have the following cases, see Figure]^ left side. 

1. Pi is not matched. This contributes Ck matchings. 

2. Pi is matched to p 2 i with 1 < i < A:. This contributes Yli<i<k Ci-ibk-i matchings. 

3. Pi is matched to P 2 i-i-i with I < i < k, p 2 i and p 2 i-i -2 are not matched to each other. This 
contributes Yli<i<k matchings. 

4. Pi is matched to P 2 i-i-i with 1 < i < A; — 1, p 2 i and p 2 i +2 are matched to each other. This 
contributes X^i<j<fc_i Ci-ibk-i-i matchings. 

This yields 

bk — Cfc -|- ^ ^ Ci—ibk—i T ^ ^ Oii—lCk—i T ^ ^ Cj—l&fc—i—!• (4:) 
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Recursion for Cfc. Clearly, cq = 1. For A: > 1 we have the following cases, see Figureright 
side. 

1. Pi is not matched. This contributes ak-i matchings. 

2. Pi is matched to p 2 i with 1 < i < k. This contributes Yli<i<k Ci-iCk-i matchings. 

3. Pi is matched to P 2 i+i with 1 < i < fc — 1, p 2 i and p 2 i +2 are not matched to each other. 

This contributes ai-ia^-i-i matchings. 

4. Pi is matched to P 2 i+i with 1 < i < k — 1, p 2 i and p 2 i +2 are matched to each other. This 

contributes Ci-iCk-i-i matchings. 

This gives 


Ck — afc-i + 7 ^ Ci-iCk-i + 


^ ^ ai—iQ-fc—i—1 T 


E 


Cl—lCfc_2—1. 


After simplifying equations ([3||^, we obtain; 

fc-l fc-l k-2 k-2 k-3 

0>k — Cfc Ck—i + ^ ^ biCk—i—i + ^ ^ QO/c—1—i T 2 ^ ^ 6jCfc_2— i T ^ ^ CiQk—2—i T ^ ^ 

2 = 0 2=0 2=0 2=0 2 = 0 

k—1 k—1 k—2 

bk — Cfc + ^ ^ 1—1 “h ^ ^ ^{C-k—l—i H“ ^ ^ 2— 2 

2 = 0 2 = 0 2=0 

fc-l fc-2 k-2 

Cfc C^fc—1 ^ ^ ^ ^ ^ ^ ^i^k—2—i 

2=0 2=0 2=0 


(5) 
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Figure 6: The cases in the recursions for bk and c^. 


We translate these equations into generating functions and obtain the following system, where 
we write A, B, C for A(x), B{x), C{x): 

A = C{{1 — x) x(l + x)A + x(l + x)‘^B) 

B = C(1 + xA + x(l + x)B) 

(7 = 1 + xA + x^A^ + x{l + x)C‘^ 

We eliminate A and B from this system and find that C satisfies the equation 

1 — (1 + 3x + 5x‘^)C + x(5 + 8x + + 9x^)(7^ — 8x^(1 + x)(l + x + x^)C^ + 

+ 4x^(1 + X + x^)(l + x)^(7‘^ = 0, (6) 


and that A and B are related to C as follows: 

(7(1 - X + 2x2(7 + 2x3(7) 


A = 


B = 


1 - 2x(7 - 2x2(7 
(7(1 - 2x2(7) 


1 - 2xC - 2x2(7 

Equation ® has four solutions. Only one of them can be written as a formal power series: 


c 


2(1+X+x3) —.^2(1+X+x3) ^1 —2x —8x 2—3x3 + (l+x)-y/ 

4x(l+x)(l+x+x3) 


The other three solutions have different combinations of signs before the two square roots. For 
those combinations, the numerator has a non-zero constant term, and this cannot balance the 
absence of a constant term in the denominator. For the series (7(x) given above, the singularity 
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closest to 0 occurs in /i = |, one of the roots of 1 — 9x — 3x^. It is a square root singularity, 

and there is no other singularity with the same absolute value; thus, by the exponential growth 
formula 0 Thm. IV.7] and a transfer theorem jSl Thm. VI.1], the asymptotics of the sequence 
is Ck = 0((l/^)''A;-3/2) with 1/// = (\/93 + 9)/2 ss 9.3218. 

Since Ck counts matchings of 2k points, it follows that the number of down-free matchings of 

SZZC„ of this kind is 0(A”/n^/^), where A = y^l/^ = 'Y/(\/ii”-|-~9)/2 « 3.0532. It is easy to see 
that the same bound holds for all kinds of zigzag chains: for the proof, note that a zigzag chain 
of kind C with 2k points includes a zigzag chain of kind A with 2A: — 1 points and is included in 
a zigzag chain of kind A with 2A: + 1 points; similarly for kind B. 

Finally, it follows from Theoremj^that the number of perfect matchings of DZZCn (of either 
kind) is fI(A"’/n^) and 0(A"'/n^). This proves Theorem]^ 

5 r-chains without corners 

5.1 Definition of r-chains with and withont corners 

In the following two sections we deal with further generalizations of the double chain. An upward 
single chain will be called an arc. As usual, the size of an arc is the number of its points. Recall 
that three points with distinct x-coordinates are in upward position if they form a clockwise 
oriented triangle when sorted by x-coordinate. 

We dehne an r-chain (with corners) with k arcs, to be denoted by CH(r, fc), see Figure j^a) 
for an example. It consists of k arcs of size r + I, the rightmost point of the ith. arc (I < i < k — \) 
coinciding with the leftmost point of the [i + I)st arc, so that any three points are in upward 
position if and only if they belong to the same arc. An r-chain CH(r, k) has rk + 1 points. As 
a special case, a simple (downward) chain is a I-chain, and an even zigzag chain of odd size is 
a 2-chain. 

One can construct an r-chain CH(r, k) with k arcs as follows: 

• Take k + 1 points Vq, Vi, V 2 ,..., 14, sorted by x-coordinate, in downward position. These 
points will be called the corners. 

• For each i = 1,2,... ,k, add r — I points on the segment Vi-iVi. 

• Replace each segment Vi-iVi by a very flat upward circular arc through Vj_i and Vi. Move 
the r — I points from the segment vertically upwards so that they lie on this circular arc. 
The radius of the circular arc must be sufficiently big so that the orientation of triples of 
points that do not lie on the same segment is not changed. 

We shall often use a compact schematic drawing of r-chains as in Figure [^b). In such drawings 
we have to draw some matching edges as curved lines rather than as straight-line segments. 

The class of (double) r-chains was earlier used for finding lower bounds on the maximal 
number of triangulations (tr) of point sets in the plane. Garcia, Noy, and Tejel [B] showed 
that tr(DCn) = 0*(8”). Aichholzer, Hackl, Huemer, Hurtado, Krasser, and Vogtenhuber [T] 
improved this bound by showing that tr(DZZCn) = 0*(8.48”). This result was further improved 
by Dumitrescu, Schulz, Sheffer, and Toth [3], who showed that a double 4-chain of size n, denoted 
in their work by D{n,3^^^), has H(8.65"') triangulations. 

Now we dehne a variation of this structure whose analysis is easier. An r-chain without 
corners with k arcs, denoted by CH*(r,/c), is a set obtained from CH(r -|- l,k) by deleting the 
corners. It consists of rk points. See Figure for an example. In this section, we will analyze 
r-chains without corners, and we will hnd precise asymptotic estimates for the number of down- 
free matchings. In the next section, we will turn to r-chains with corners. They give even 
stronger lower bounds, but the analysis will not be so precise. 
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Figure 7: A 5-chain (with corners) with six arcs: (a) a precise drawing; (b) a schematic drawing. 


5.2 Matchings with runners 


Consider a matching of X = CH*(r, A:). We want to build down-free matchings incrementally 
from left to right by adding one arc at a time. If we cut such a matching between two arcs, 
then we possibly have some edges cut into two “half-edges”, which we call runners. (In botany, 
runners are shoots that connect individual plants.) A runner can be formally defined as a marked 
vertex. Such a vertex must not be matched by “proper” edges and must be visible from above. 
These requirements ensure that, in the course of the incremental construction, two runners can 
be joined into one edge. Runners are visualized as half-edges that have one endpoint in X and 
the other end dangling, see Figure [^a). Note that it is not assigned in advance whether a runner 
will be matched to the left or to the right. 

A matching which possibly has runners will be called a p-matching. Extending our previous 
definition of free points, we call a point free in a /o-matching if it is neither matched by a “proper” 
edge nor marked as an endpoint of a runner. A p-matching is down-free if all free vertices are 
visible from below. 

In the course of the recursive construction of down-free p-matchings, runners from different 
arcs can be matched as follows. Let A and B be two r-chains without corners, and let Ma 
and Mb be down-free p-matchings of these sets. We place B to the right of A. If Ma has j 
runners and Mb has /3 runners, then for each I in the range 0 < £ < min{j,/?} we can match, 
in a unique way, the rightmost I runners of Ma to the leftmost i runners oi Mb. The obtained 
p-matching M is also down-free; the runners which were not matched in this procedure remain 
runners in M ; the number of such runners is j -|- /3 — 2i. Conversely, each down-free p-matching 
of AuR can be obtained by this procedure from two uniquely determined down-free p-matchings 
of A and B. Figure |^b) shows an example with j = 4, /3 = 3, ^ = 2. 

We summarize these observations for the special case that we will use in the recursive con¬ 


struction of r-chains: adding one new arc to the right of a given r-chain, see Figure 10 


Observation 9. Let X = CH*(r, k) be an r-chain without corners with k > 1 arcs. Let B be the 
rightmost arc of X, and let A = X\B. Let Ma be a down-free p-matching of A with j runners, 
and let Mb be a down-free p-matching of B with /3 runners. For each 0 < £ < min{j,/3} there 
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(a) 

(b) 

• • •• •• •• •• •• 

• •• •• •• •• •• • 

Figure 8: A 4-chain without corners with six arcs: (a) precise drawing; (b) schematic drawing. 





Figure 9: (a) A down-free p-matching Ma with four runners of A = CH*(6,3). (b) Combining 
Ma with a down-free p-matching with three runners of B = CH*(6,1). 


exists a unique down-free p-matching Mx,e of X obtained by matching the rightmost i runners 
of M with the leftmost i runners of N. The number of runners in Mx/ is i = j (3 — 2L 

Conversely, each down-free p-matching M of X can be obtained in this way from uniquely 
determined matchings Ma and Mb (of A and B) as above. If M has i runners, Ma has j 
runners, and Mb has /3 runners, then the number of edges obtained by matching of pairs of 
runners is (. = {j -\- (3 — i)/2. 

For k = 1, this observation holds trivially: A is empty, and the only possibility is j = I = 
0, (3 = i. From the above relations between the parameters i,j,l3,i, one can work out the 
constraints on the possible values of /3 for given i and j: The equation i = j -\- f3 — 2i together 
with 0 < i < min{j, (3} implies that (3 must satisfy \i — j\ < (3 < i + j and (3 = i — j (mod 2). 

5.3 Recursion for matchings with runners in r-chains without corners 

Denote the number of down-free p-matchings of CH*(r, k) with i runners by z^{r) or simply by 
z^, since we will regard r as fixed. Obviously, the down-free matchings of A = CH*(r,/c) are 
the down-free p-matchings without runners. Since the size of X is rk, the growth rate for the 
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j * 



Figure 10: Runners in a recursive construction of a /O-matching of an r-chain without corners. 


number of its down-free matchings is ^ 

For A: = 0 we have Zg = 1 and = 0 for z > 0. The numbers zj for a single arc will serve 
as a basis of the recursion. They are determined in the following proposition. 

Proposition 10. 1. The number of down-free matchings (without runners) of a single arc 

of size r is 

2. The number of down-free p-matchings of a single arc of size r that have i runners is 


4 = 44 ) = Q 


r — I 

L(r - i)/2\ 


i, [(r-z)/ 2 j, [(r-i)/ 2 ] 


Proof. 1. For the first equation, let f{x) = 44)^^ be the generating function for the 

number of such matchings in terms of the size r of an arc. We will show that f{x) satisfies the 
equation 

fix) = • c(x^) • fix) + l) , (7) 

where c(x) = (1 — Vl ~ 4x) /2x is the generating function of the Catalan numbers. Therefore, 
we have 

f4) = 


1 


1 — X — x^cix"^) ’ 

and this is known to be the generating function for (Lr/2j) ig A001405]. 

To see why 0 holds, consider the leftmost matched point p (if there is any). Suppose that 
p is matched with q, see Figure pT] for illustration. Then all points to the left of p are free, which 
contributes 1/(1 — x) to the generating function. The points between p and q are not visible from 
below and, therefore, they are matched by a perfect matching; this contributes c(x^). Finally, 
the points to the right of q are matched by a down-free matching, whose generating function 
is again fix). The factor x^ accounts for the two points p and q, and the additive term -|-1 
accounts for the case that p does not exist. 

We give another proof - a bijective one. For a given matching, we mark the left and right 
endpoints of each edge by L and R, respectively. We leave the free points unmarked for the 
moment. Then the non-crossing matching can be reconstructed from the labels: We traverse 
the points from left to right, and we match each R that we meet with the closest previous 
unmatched L. Moreover, since the matching is down-free, free vertices can only appear when 
there are no previous unmatched L-vertices. Now we label the free points: If there are 7 free 
points, we label the first [ 7 J free points by R and the last [ 7 ] free points by L, see Figure 12 
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o 




Figure 11: The leftmost edge pq in 


the proof of Proposition 10|1 


for illustration. The free points marked R can be recovered in a left-to-right sweep as those 
R-vertices for which we find no previous matching L-vertex in the above procedure. The free 
points marked L can be recovered similarly in a right-to-left sweep, and finally, the matching 
among the non-free points can be found as described above. Thus we have established a bijection 
with sequences of length r over the alphabet {L,R} that contain [r/2j R’s. 


• • 

• • 

RRLLL 



LLRRL 


/ •. 

LRRLL 

/ 

LRLRL 



LLLRR 


• • 

RLLRL 

\ 

LLRLR 


; \ 

RLLLR 

/ \ 

LRLLR 


Figure 12: The coding of down-free matchings in the second proof of Proposition 10|1 


2 . Let us turn to the second equation. Once we choose i endpoints of runners, the whole 
matching is down-free if and only if its restriction on the remaining r — i points is down-free. 
Therefore, the result follows directly from the hrst part. □ 

Now we find a recursion for zf,k>l. 


Proposition 11. For fixed r, we have the recursion 



E 


CLijZj 


A-1 


j>0 


( 8 ) 


with coefficients 


°'ij - + ^\i-j\+2 ^ ^mm{r*,i+j}^ ( 9 ) 

0</3<r, 

li-jl<h<i+j, 

^=i—j (mod 2) 

where r* is r or r — 1 and has the same parity as i — j. 

Proof. For /c = 1, Q can be verified directly. Assume now X = CH* (r, k) with /c > 1, let i? be 
the rightmost arc of X, and let A = X \ B. For each j > 0 and each possible /3 we will find the 
number of p-matchings of X with i runners whose restriction to A has j runners and restriction 
to B has (d runners. By Observation p-matchings of A and B and the values of i, j and /3 
determine uniquely an p-matching of X. Therefore p-matchings of A and B with (respectively) 
j and /? runners contribute ■ z^ p-matchings of X with i runners. 

For given i and j, the bounds \i — j\ < fd < i + j and the restriction fd = i — j (mod 2) given 
in are explained in the remark after Observation □ 
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5.4 Analysis of the recursion 

For each k > 0, denote Vk = {zq, , -zf, • ■ • In particular we have vq = (1,0,0,0,... 

Consider the inhnite coefficient matrix A = (aij)jjgNp with given by ([^. By Proposition [IH 
we have Av^-i = for each fc > 1. One easily verifies that the matrix A has the following 
properties: 

• A is symmetric. 

• A is a band matrix of bandwidth r: for \i — j\ > r we have aij = 0 . 

• The entries of the first row and column are any = ooi = zj = (0 ([(r- 7 v 2 j)' 

• For f + j > r* we have aj+ij+i = Uij. That is, the diagonals - sets of entries with hxed 
q '■= j — i, \q\ < r - stabilize starting from the entry a(r*-g)/ 2 ,(r*+<j)/ 2 - For these entries we 
have: 

aij — ttijy-q — ^ ^ zp. (fO) 

0=q (mod 2) 

In particular, starting from the rth row (resp. column), the rows (resp. columns) are shifts 
of each other, and therefore, have the same sum of elements. 


The elements in the upper-left corner {i+j < r*) are positive and smaller than the elements 
in the same diagonal after stabilization - since in this case we have a partial sum of ( 10 ). 


For example, for r = 5, the matrix is 


no 

30 

30 

20 

5 

1 

0 

0 

0 

0 

0 

...X 

30 

40 

50 

35 

21 

5 

1 

0 

0 

0 

0 


30 

50 

45 

51 

35 

21 

5 

1 

0 

0 

0 


20 

35 

51 

45 

51 

35 

21 

5 

1 

0 

0 


5 

21 

35 

51 

45 

51 

35 

21 

5 

1 

0 


1 

5 

21 

35 

51 

45 

51 

35 

21 

5 

1 


0 

1 

5 

21 

35 

51 

45 

51 

35 

21 

5 


0 

0 

1 

5 

21 

35 

51 

45 

51 

35 

21 


0 

0 

0 

1 

5 

21 

35 

51 

45 

51 

35 


0 

0 

0 

0 

1 

5 

21 

35 

51 

45 

51 


0 

0 

0 

0 

0 

1 

5 

21 

35 

51 

45 


i : 











•• 


( 11 ) 


The column sum after stabilization of the columns, that is, starting from the (r -|- l)st 
column, is as follows: 


Ay. 


^(i + l)4 = ^(i + 1) 
2=0 2=0 



r — i \ 


E 

i=0 


(^+ 1 ) 


f, L(r-f)/ 2 j,r(r-f)/ 2 l 


( 12 ) 


In the spirit of the Perron-Frobenius theorem for non-negative stochastic matrices, one can 
expect that A^ is the growth rate for {zQ)k>o- We will prove the this is indeed the case by using a 
result by Banderier and Flajolet [3] about enumeration of certain kinds of colored lattice paths. 

Proposition 12. For fixed r, we have Zq = Q*{{Xr)^). 

Note that the superscript k in the left-hand side denotes an index, whereas in the right-hand 
side it is a power. 
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Proof. We begin with some notion for lattice paths. Families of lattice paths are usually defined 
by indicating a starting point - normally ( 0 , 0 ) ~ and a set of possible moves of the form (l,/3). 
For many familiar families it is additionally required that the paths never go below the x-axis 
and/or end at the x-axis. The paths that start at (0,0) and satisfy both these restrictions are 
called excursions. For example, Motzkin paths [HI A001006] are excursions that use the moves 
( 1 , 1 ), ( 1 , 0 ), ( 1 ,- 1 ). 

In a more general setting, a set of possible moves may depend on the point reached by a 
path. Moreover, each move (l,/3i) starting in certain point can have a non-negative integer 
multiplicity mi. This is sometimes expressed by saying that these are copies of the same move 
that are distinguished by mi different “colors”. 

In summary, to each lattice point (a, h) we assign a rule — a set of moves that can be used 
for the next step once a path reached this point, together with multiplicities. It is assumed that 
for each lattice point the number of moves with non-zero multiplicity is finite. Note that one 
can express the condition of non-crossing the x-axis in terms of such rules: one has to require 
that for each point (a, b) there are only moves (1, /3) with /3 > —h. 

Consider now the case that all points that lie on the same horizontal line have the same rule. 
Namely, for y = j and i > 0 we denote by dij the multiplicity of the move (1, i —j) at (any) point 
{a,j). We collect these data in the infinite matrix D = (djj)ijgNo- Let u = (1, 0, 0,...)"'". Then 
the number of paths that start at ( 0 , 0 ), do not cross the x-axis, and end at a point (a, b) is equal 
to the 6 th component of - this follows directly from matrix multiplication. In particular, 
the upper-left entry of is the number of excursions of length k, which we will denote by 
Px{D,k). The quantity in which we are interested, the number Zq of down-free matchings, is 
then given by Zq = Ex(A, k), where A is the coefficient matrix given above (10). 

Suppose now that we have an even more restricted case: all points have the same rules; yet 
still we want to consider only paths that remain weakly above the x-axis, so we exclude the 
moves that violate this requirement. For such families, a result of |3l Theorem 3] can be applied. 
It states that the number of excursions of length k with moves {(l, 6 i),(l, 62 ),...,(l, 6 m)} and 
associated multiplicities rci,..., Wm, is of the form Q{C^/k^^‘^), where the base C of the expo¬ 
nential growth is determined as follows: For the Laurent polynomial P{u) = WjU^y let r 

be the unique positive number such that P'{t) = 0; then C = P{t). The situation is particularly 
easy for families with a symmetric set of moves, that is, if ( 1 , 6 ) is a move then ( 1 , — 6 ) is also 
a move with the same multiplicity, or equivalently, P{u) = P{u~^). In this case, r = 1, and 
consequently, C = P{t) = Yl'JLi Wj. 

The situation for our matrix A is very similar to this case, except that the first r — 1 
horizontal lines of the lattice follow different rules, in accordance with the fact that the first 
r — 1 rows of A are different from the others. However, this does not affect the asymptotic 
growth rate. Indeed, let us look at the matrix A' in which the first r rows and columns of 
A have been removed. It coincides with A for i -|- j > r, but the rule aj+ij+i = Oij holds 
for all entries - also in the upper-left corner. Since A < A' element-wise, we clearly have 
Ex(A, k) < Ex(A', k) = 0(A^/fe^/^). To see that we have a lower bound of the same asymptotic 
form, consider only those excursions that start with the move ( 1 , -|-r), end with the move ( 1 , —r), 
and never go below level r. The intermediate part of the excursion is governed by the matrix A 
from which the first r rows and columns have been removed, which coincides with the matrix A'. 
Thus Ex(A, k) > Ex(A', k-2) = Q{X^/k^/^). □ 


5.5 Asymptotic growth constants 


Since A = CH*(r, k) has n = rk points, it follows from Proposition 12 that the growth rate for 
the number of down-free matchings of the r-chain without corners of size n is \/A/. In order to 


estimate A,., we note that the expression (12), when the factor (i-|-l) is ignored, counts partitions 
of r elements into three subsets (the latter two being of almost equal size). The total number 
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of such partitions is 3^. Hence, < (r + 1)3^, and \/A7 converges to 3. Computations suggest 
that the maximum of \/A^ is obtained for r = 11: ^/An = V240054 ~ 3.0840; after that it 
apparently decreases monotonically to 3, see the left part of Table for the first few values. To 
prove that r = 11 gives indeed the maximum, one estimates that < 3^r + 1 < 3.0838 for 
r > 191, and the finitely many values up to r = 190 can be checked individually. This completes 
the proof of Theorem 


In order to find a more precise estimate for A^-, we notice that the middle expression in (12) 
expresses \r as the binomial convolution of the sequence of natural numbers and the sequence 
,). It follows that the exponential generating function for {Xr)r>o is 


' m 


(1 + x) {Io{2x) + Ii(2x)), 


where Io{x) and Ii{x) are the modified Bessel functions of the first kind. From this we can 
conclude that the sequence (Ar)r>o is the sum of the sequence A005773 and a shifted copy of 
A132894 in [8]. The ordinary generating function for this sequence is then 

1 / 1 — 2x — x^ 

2x V(1 + xYG{l — 3x)3/2 

and it follows from the exponential growth formula that A^ = 0(3'’r^/^). By Theorem this 
number is also the growth rate of the number of perfect matchings for the corresponding double 
structure. 



6 r-chains with corners 

6.1 Recursion 

In this section, we will treat r-chains with corners, but we will simply refer to them as r-chains. 
The analysis of these r-chains is more complicated due to the fact that the corners belong to two 
arcs. As before, we will incrementally build the r-chain and estimate the number of matchings 
of the r-chain with k arcs, which possibly have runners. We extend the notions of runners, free 
points, and p-matchings to r-chains with corners in the obvious way. 

We cut a down-free p-matching M of CH(r, k) to the right of I4_i - the rightmost point of 
the {k — l)st arc - and obtain two down-free p-matchings: the first, Ma, of A - the set consisting 
of the first k — 1 arcs of CH(r, k)] and the second, Mb, of H - the rightmost arc of CH(r, k) 
without the point I4_i. See examples in Figures 13-18 Note that in the case of r-chains with 
corners a runner incident to I4_i, upon adding B on the right, can be also connected to a point 
of H: in such a case we say that it is matched internally. 

We distinguish whether M has a runner incident to Vk or not. Let be the number of 
down-free p-matchings of CH(r, k), where 14 has a runner and in addition there are i runners. 
Let Ph be the number of down-free p-matchings of CH(r, k), where 14 has no runner and there 
are i runners. (C stands for “corner”, F for “free”.) For A: = 0, there is a single vertex, and we 
have Cg = Tq = 1 and = 0 for all i > 0. The number that we are interested in, the 

number of matchings in CH(r, k), is Fq. Next we find recursive expressions for and for 


Recursion for Cf. For Cf, the new corner 14 has a runner and is not available for receiving 
edges from the left. Thus for the formulae below, it can be treated as if it were not present in 
the fcth arc. We have the following three cases: 


1. (Figure 13 ) The previous corner 14_i has a runner which is not matched internally in 


the A:th arc. Suppose there are a runners originating in the fc-th arc, in addition to that 
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without corners 

with corners 

r 



Mr 

Tr 

1 

3 

3 

ill) 

3 

2 

9 

3 

(7 i) 

3.0532 

3 

28 

3.0366 

(^?J) 

3.0711 

4 

87 

3.0541 

l' 31 28'i 

Ue 59 1 

3.0819 

5 

271 

3.0662 

to 

OCD 

00 00 

3.0877 

6 

843 

3.0735 

( 301 2711 

1 632 572 ) 

3.0909 

7 

2619 

3.0783 

( 933 8431 

1 1952 1776 ) 

3.0925 

8 

8123 

3.0812 

( 2885 2619 1 

1 6022 5504 ) 

3.0930 

9 

25153 

3.0829 

( 8907 81231 

1 18550 170401 

3.0929 

10 

77763 

3.0837 

( 27457 251531 
i 57071 526101 

3.0923 

11 

240054 

3.0840 

( 84528 777631 
i 175381 1622911 

3.0915 

12 

740017 

3.0839 

( 259909 240054 1 

1 538386 499963 1 

3.0904 

13 

2278329 

3.0835 

( 798295 7400171 

1 1651140 15383121 

3.0893 

14 

7006093 

3.0829 

( 2449435 2278329 1 
i 5059251 4727764 1 

3.0880 

15 

21520872 

3.0822 

( 7508686 70060931 

1 15489221 14514779 1 

3.0867 

16 

66039651 

3.0813 

( 22997907 21520872 1 
f 47384904 44518779 1 

3.0854 

17 

202462113 

3.0804 

( 70382811 660396511 
t 144857454 136422462 1 

3.0841 

18 

620164491 

3.0794 

( 215240265 202462113 1 

1 442540653 417702378 1 

3.0828 

19 

1898109900 

3.0785 

( 657780918 6201644911 

1 1351126551 1277945409 ) 

3.0815 

20 

5805127269 

3.0774 

( 2008907469 1898109900 1 

1 4122747150 3907017369 J 

3.0803 


Table 1: Summary of results for r-chains without and with corners, for 1 < r < 20. For r-chains 
without corners, Xr is the row sum of the matrix A (Section |5.4[ ), and \/^ is the growth rate for 
pm. For r-chains with corners, the condensed coefficient matrix Mr is derived from the recursion 


(Section 6.4), and Tr, the r-th root of its dominant eigenvalue, is the growth rate for pm. In 
both cases, the values for r = 1 and r = 2 reproduce the known bounds. Indeed, a 1- and a 
2-chain withont corners, as well as a 1-chain with corners, is just a downward chain, and thus 
the growth rate of 3 agrees with Theorem A 2-chain with corners is a zigzag chain, and thus 
T 2 ~ 3.0532 agrees with Theorem 


^k—1 


(13) 


originating in V^-i- These runners can be only matched to the right. The contribution to 
C‘is 

E 

0<Q:<min{r— 1 , 2 —1} 

where 

G — l\ / r — 1 — a 
a )\f{r - I- a)/2\. 

The expression for is similar to from Proposition 10|2 but here we have only r — 1 
points: all the points of the A:th arc, excluding the corners. 

2. (Fignre [14} ) I4_i has no runner. This possibility contributes 


ZcL — 


E E 

j>0 \i—j\<<y-<i-\-j 

OL=i—j (mod 2) 
0<a<r—1 


k-l 


(14) 
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the fcth arc 

Figure 13: Case 1 in the recursion for C^: V^-i has a runner not matched internally in the kth 
arc. 

This formula (as well as some of the formulae in the following cases) has the same pattern 
as Q, with appropriate changes. 



the fcth arc 

Figure 14: Case 2 in the recursion for C^: 14_i has no runner. 


3. 


(Figure 15 
is 


14_i has a runner matched internally in the k-th arc. The contribution to 

E E (15) 

j>0 

a=i—j (mod 2) 

0<Q(<r—1 


where 


la 



f r — a \ 
[[{r-a)/2\) 


r — 1 — a \ 
L(r-l-a)/2jj_ 


/r — 1\/ r — 1 — a \ 
V a )[[{r-2-a)/2\)- 


In the expression for la, the first factor counts the choices of a runners from the r — 1 
points. In the second factor, we subtract from all down-free p-matchings on the remaining 
r — a points (including 14 _i) those in which Vk-i is unmatched, which is the same as 
down-free p matchings on r — 1 — a points. 
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the fcth arc 

Figure 15: Case 3 in the recursion for C^: 14_i has a runner matched internally in the /cth arc. 


Cf is the sum of the three expressions (13 15). 


Recursion for F^. For Fj^, we have again three cases: 

1. (Figure [T^) 14_i has a runner not matched internally in the fcth arc. In this case, all the 
additional a runners originating in the interior of the k-th arc must be matched to the 
right. Vfc is either free or matched internally to the left. The contribution to Fj^ is 

E (16) 

0<Q;<min{r— 1 , 2 —1} 


where 


Wry = 


r — 1 
a 


r — a 

L(r - a)/2j 


is again similar to from Proposition 10|2 but here we have r 
because no runner originates from 14 . 


1 in the first factor 



the fcth arc 

Figure 16: Case 1 in the recursion for Fj^: 14 -1 has a runner not matched internally in the /cth 


arc. 
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2. (Figure [I t}) 14 has a runner connected to a point of A\{14_i}. In this case, all a runners 
originating in the /c-th arc must be matched to the left. The contribution to is 


E + z«d+-i+o)- (17) 

0<Q<r—1 


The two terms - with ^ and with ^ - correspond to the subcases where I4_i is 
internally matched or, respectively, not matched to a point of the fcth arc. 



the kth. arc 


Figure 17: Case 2 in the recursion for Fj^: 14 is connected to a point to the left of I4_i. 


3. (Figure [T^) I4_i has no runner, and 14 has no runner matched to a point of A \ {I4_i}. 
The contribution to Fj^ is 


where 



E E iUo.C^-^ + W^F^-^), (18) 

j>0 

a=j—i (mod 2) 

0<Q<r—1 


/ r + l — a \ / r — a \ 

[[{r + l-a)/2\)~[[{r-a)/2\)_ 



f r — a \ 
U(r-l-a)/2jJ 


The two terms correspond to the same possibilities as in the previous case. The factor Ua 
is similar to Iq in the third case for Cf, but here we count the down free p-matchings of 
the whole kth. arc with both its corners, hence we have r + l instead of r in the second 
factor. 


Fk is the sum of the three expressions (16 18). 


6.2 Analysis of the recursion 

The expressions above imply a coupled mutual recurrence between two sequences of vectors = 
{C^, C^, .. .y and = {F^F^ . The initial values are = F^ = (1, 0, 0,.. .)^. 

and F^ are expressed in terms of and F^~^ as follows. For i > r, we have: 


+ = E 

l3=—r 


CC(jk-l 


i+P 

0=—r 


^CF pfc-1 


+ = E 


a/3 '-'i+13 + 


E ^FF pfc-1 
^i+l3 > 

l3=-r 


(19) 
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the fcth arc 

Figure 18: Case 3 in the recursion for F^: V^-i has no runner and 14 is not connected to a 
point to the left of 14 _i. 


where the numbers , aP^■, are to be read out from the expressions in Section 6.1 

For the small indices i < r, we have irregularities, like for r-chains without corners: The co¬ 
efficients in (19) must be replaced by smaller coefficients which depend also on h In matrix 


notation, the recursion is written as 


_ j^CC ^k—1 _j_ j^CF jpk—1 


( 20 ) 


where there are four band matrices , A^^, A^^, A^^ of bandwidth r, similar to the matrix 
A from ( [IT| ). 



Figure 19: The recursion (19) gives the number of paths on this network. The neighborhood of 


a typical vertex Ci is shown in a schematic way. 


This system can be interpreted as a set of lattice paths on a two-layer lattice, see Figure IT 
We have a row of nodes Cq, Ci, 6 * 2 ,... and another row of nodes Fq, Fi, F 2 ,... immediately below 
it. The possible jumps and their multiplicity depend only on the row, with irregularities close to 
the left edge. (In this representation, the lattice paths considered in the proof of Proposition [T^ 
in Section [5.3| correspond to walks on a ray 0,1, 2,.... The x-coordinate of the two-dimensional 
lattice in Section [5.3| is now represented as time. 

We are not able to provide as precise estimates for the growth rate as for chains without 
corners, where we had a single recursion. One would expect a similar behaviour. However, we 
can still pin down the base of the exponential growth as an eigenvalue of an associated 2 x 2 
matrix. 

First, we can get rid of the irregularities but cutting of the first r rows and columns of 
the coefficient matrices. As for the case of a single matrix, this does not affect the asymptotic 
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growth. We can now assume that the diagonals are constant, and the recursion (19) holds for 
all i, with the convention that and Pj~^ in the right-hand side are taken as 0 for j < 0. 

For better readability, we will now replace the vectors and by more generic names 

and 


Vi = 


r 

E 

P=—r 

r 

E 

[3=—r 


XX fc-l , 
0/3 + 


r 

E 

j3=—r 




( 21 ) 


VX fc-1 , 
®/3 ®j-r/3 + 


0/3 y/^/3, 

l3=-r 


for all i, with the understanding that quantities and with negative subscripts j on 
the right-hand side are regarded as zero. 

Our analysis below is valid under the condition that the coefficients of this recursion don’t 
exhibit a tendency to favor larger or smaller indices, or in other words, that the Markov chain 
associated to the system does not systematically drift to the left or to the right. (In the one- 
vector recursion analyzed in the proof of Proposition this no-drift condition was not an issue 
because the set of moves was symmetric.) To formulate this condition precisely, we have to set 
up some notation and establish some terms. 


Let us denote the coefficient sums in the terms of the recursion (21) as follows: 

r 

1 '-^= j : 




r 

E' 

/3=—r 


.XX 


= 


r 

E 

l3=—r 


®/3 ) 


0-13 y 


X'-' = 


E 


YY 


l3=—r l3=—r 

These numbers are the column sums of the coefficient matrices after stabilization. These sums 
form the condensed coefficient matrix 

^^XX ^XY'' 

^YX JYY 


( 22 ) 


Let M denote its dominant eigenvalue. Let {px,Py) be the corresponding left eigenvector and 
(ttx, vTy)"'^ be the corresponding right eigenvector, with the normalization px+PY = xx+xy = 1. 
Since the matrix is positive, these two vectors are positive. 

We define the total group-to-group jump sizes of the system: 




E 


„xx 

0/3 


13, D-3'- = 


E' 

p=—r 


,XY 


13, D''^ = 


E' 

/3=—r 


,YX 


13, O’"'' = 


E' 

l3=-r 


,yy 

A 




The weighted total jump size D of the system is then defined as follows: 

D = pxxxD^^ + pxt^yD^^ + pyttxD^^ + py^yD^^ 


(23) 


^ (« .) --) (-) 


Now we can state the main result of the analysis. 

Theorem 13. Suppose the system (21) has non-negative coefficients, and the weighted total jump 
size D is zero. Assume that the coefficients , 0 ^^, 0 ^^ are positive for /3 = —1,0,1. 

Let M be the dominant eigenvalue of the condensed coefficient matrix (22). Then 

xg = 0(M"), = 0(M"), 


and 

xt = n{{M-ef), y'^ = n{{M-ef) 

for every e > 0. 

Since the proof is quite substantial, we devote a separate section to it. 
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6.3 Proof of the theorem about mutually coupled recursions 

We will transform the problem to a recursion in which the left eigenvector is {px, Py) = (Ij 1); 
and thus the column sums of the coefficient matrix (after stabilization) are constant. We achieve 
this by rescaling the vectors x and y to = pxx^ and = PyUi ■ Clearly, the asymptotic 
growth of X and y is unaffected by this multiplication with a constant. For these new vectors, the 
coefficients of the recursion change to = px/ Py- and = py/px- ■, while = 
are unchanged. Consequently, the first column sum of the condensed coefficient 


matrix ([22|) becomes + Py/px ■ A^^ = {px ■ A^^ + py ■ A^^)/px = {Mpx)/px = M, 


and similarly for the second column. Theorem |13| follows therefore from the following theorem, 


which is a special case of Theorem 13 with the additional assumption that the matrix (22) has 
constant column sums. 


Theorem 14. Suppose the system (21) has non-negative coefficients and constant column sums 


M = A^^ +A^^ =A^^ + A 


YX 


TXY 


iYY 


Suppose that 


vrx + vry = 0, 


(24) 


(25) 


where (ttxj'^y) is a right eigenvector of the matrix ( [2^ with eigenvalue M. Suppose further 
that the coefficients , a^^, are positive for j3 = —1,0,1. Then 


and 


for every e > 0. 


x^ = O(M^), y^ = O(M^), 


x^ = Q((M-e)^), y^ = Q((M-e)^) 


Theorem 14 is formulated in terms of the original recursion (21), but it must be applied to x 


and y instead of x and y in order to prove Theorem 13 Our rescaling modifies group-to-group 
jump sizes in the same way as the coefficients: = px! Py ■ , etc.; the eigenvectors of 

the modified condensed coefficient matrix are (ffxjTry) = (px^ry, pyvry) and {px,Py) = ( 1 ) 1 ) 
(without normalization), and with these substitutions, the condition that D from ( |23[ ) is zero 
translates into (|25|), after erasing the tildes. This concludes the proof of Theorem fTsI □ 


Proof of Theorem\14\ The upper bound is easy: by summing all equations of (21), one sees 
that X)i>o hi grow at most by the factor M in each iteration, since the column 

sums of the coefficient matrix are bounded by M. It follows that XQ,yQ < + YliVi ^ 

Let us now turn to the lower bound: To have a compact notation for the linear operator 


expressing in the recursion (21), we denote it by <j)-. 

(x^/) = </>(x"-^ 


.k-l\ 


As an intermediate lemma, we will show that any “sub-eigenvector” with eigenvalue A is enough 
for a lower bound on the growth. 

Lemma 15. Suppose there is a pair of non-negative non-zero vectors x and y with finitely many 
non-zero elements such that the inequality 


4>{x,y) > A • {x,y) 


(26) 


holds componentwise for some A > 0. Then there is a constant K > 0 such that XQ,yQ > KX^ 
for all n G N. 
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Proof. Since (f is a monotone operator, the inequality (26) remains fulfilled if we apply 4> on the 
vectors arbitrary many times: 


> X-cj)^{x,y) 




Applying cj) on (x, y) sufficiently many times, we eventually obtain a vector whose components 
xo and yo ai'e positive, since the coefficients are positive by assumption. 

Moreover, by scaling we can obtain a vector in which these components are bigger than 1. Thus, 
we can assume that xq > 1 and yo > 1. 

Now, we find ni and K such that the following inequality holds componentwise for n = ny. 


{x\y^)>KX^-{x,y) 


(27) 


To see that this is possible, we use the assumption that are positive for 

/3 = 0 and /3 = —1. Thus, by making ni big enough, we can ensure that {x^^,y'^^) has positive 


components wherever (x,y) has positive components. We can then fulfill (27) by choosing K 
small enough. 


The inequality (27) carries over to all larger n by induction, using monotonicity of the 


operator cj) and the assumption (26). Since xq > 1 and yo > I, the desired inequalities follow 


from (27) for all n > ni. Finally, for the finitely many values n < ni, we can fulfill the inequalities 
Xo,yo > KX^ by decreasing K if necessary. □ 

Let us explain the idea for getting “sub-eigenvectors” x and y for Lemma If we wish 


to fulfill (26) for A = M, vectors x and y with constant entries will do the job. However, they 
have infinitely many non-zero entries. Thus, we aim for a smaller X = M — e, and we make 
an ansatz where the entries are determined by a concave quadratic function. This has to be 


adjusted later because the vectors have to be non-negative, and because the recursion (21) has 


some irregularities for the small values i < r. Moreover, the two coupled sequences x and y 
depend on each other in a non-symmetric way, and therefore we cannot use the same quadratic 
function for both sequences. They have to be scaled differently, and shifted horizontally relative 
to each other. 

We define the shift constant 


(5 = 


TTxD^^ TTyD 
-7ryA^>^ 


XY 


7ry(A^^ - M) 


In this definition, equality of the numerators follows from the assumption (25), which expresses 


that the weighted total jump size is zero. The denominators are equal because the column sums 


are M (24). 


We take some real parameters p and s, to be determined later, and define the quadratic 
functions hx and hy and two auxiliary vectors x and y as follows: 


hx{i) = 7rx{p- i^) 
hy{i) = xy{p - {i + S)^) 

Xi = hx{i - s) 
yi = hy{i - s) 


(28) 

(29) 


for all i € Tj. The two quadratic functions have their peaks at i = 0 and i = —5, with respective 
values pTTx and pny. These function are shifted to the right by s before they are used as 
entries of x and y. The following lemma expresses the difference of these vectors from being an 
eigenvector with eigenvalue M. 
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Lemma 16. Each of the two expressions 


Qx = M-Xi — 

Qy = M -yi- 


ap Xi+p 

,l3=-r 


^=—r 


E V X 

0/3 Xi+f) 


+ Y 


\^l3=—r p=—r 

has a constant value, independent of i, p, and s. 


(30) 

(31) 


Proof. First, we replace the quadratic function in each of the summation terms by a Taylor 
series around the weighted average point. The linear terms will then cancel, and the quadratic 
terms have a constant value. We carry this out by way of example for the sum of the terms. 
The parameters i and s always occur together in the combined term i — s, and thus we express 
our terms in terms of the parameter t := i — s. 


Y af^Xi+f}= Y aYhxiiP (d - s) = Y af^hxit + lS) 
l3=—r l3=—r l3=—r 


Let 


jjXX ^ h=-r 


E 


E 

l3=-r 


pXX 

Jxx 


denote the average jump size from group X to group X. Then we rewrite hx as a Taylor series 
in the point t + . 


hx{t + x) = hx{t + + h'x{t + D^^){x - - nxix - f 


We get 


Y af^hx{t + j3) 

^=—r 

= hx{t + Y E -^^Y f 

h h h 

= hx{t + + h'x{t + ■ 0 — const. 


We transform the other sum in the expression (30) analogously, using the average jump size 
, and we rewrite (30) as follows: 


Qx = M ■ hx{t) — hx{t + — /ly(t + + const 

= M ■ TTxip — t^) — xx{p — {t + — xy{p — {tp + 5)‘^)A^'^ + const 

= {p — t^){'nxM — ttxA^^ — vryA^^) 

+ 2t{xxD^^A^^ + ttyD^^A^^ + ttySA^^) + const 


The coefficient of {p — t^) is zero because {xx,xy) is an eigenvector, and the coefficient of t is 
zero by the definition of 5. Thus, the expression Qx has a constant value, as claimed. 
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For the expression (31), the calculation is slightly different. 

Qy = M ■ hyit) — hx{t -[- — /iy(t + + const 

= M • TTy (p — (t + (5)^) — 7rx{p — {t + A^^ — TTy (p — (t + + const 

= (p — t^)(7ryM — TTxA^^ — TTyA^^) 

+ 2t(—TTyJM + ttxD^^ A^^ + vryZi^^A^^ + TTydA^^) + const 

The coefficients of (p — and t vanish for the same reasons as above. This concludes the proof 


of the lemma. 


□ 


The vectors x and y have negative values. To get our desired vectors x and y, we will clip 
these values to 0. We determine the parameters p and s in such a way that the resulting vectors 
X and y start with a big jump from 0 to a positive value, big enough to accommodate the 
“perturbation” resulting from modifying the negative values to 0. Let e > 0 be given, and let 
K := max{Qx, Qy} be the maximum of Qx and Qy- {K is positive, but the ensuing argument 
does not depend on this fact.) We look at the sorted set of values 

I i eZ}u{(z-5)2 I i gZ} 

and find p as a positive value in this set such that the gap to the largest value which is smaller 
than p is at least K/{e min{7ri, 7r2}). Since the functions are quadratic, there must be larger and 
larger gaps as the numbers get bigger, and therefore such a value p exists. For the functions 
hx{i) = 7rx(p — and /iy(i) = 7ry(p — (i + 5)^) in (28-29), this implies that the smallest 
positive value in their range is at least K/e. Now we shift the functions horizontally such that 
positive values occur only at positive arguments, by choosing s > y/p ± <5. Finally, we clip the 
negative values and define 

Xi = max{/iy:(* — 'S),0}, xji = max{/iy(i — s),0}, 
for all z G Z. This will set Xi = iji = 0 for i < 0, in accordance with the interpretation that is 


given in (21) when these values appear on the right-hand side. 


We will show that 


(j)ix,y) > (M-e) • (x,y), 


(32) 


thus establishing condition (26) and proving the lower bound of the theorem with the help of 
Lemma [T5j 


In concrete terms, our desired relation (32) looks as follows: 


(M -e)-Xi< a^^Xi+y ^ op^Vi+y 

p=—r p=—r 

r r 

{M-s)-yi< a^Xi+y + Y 

l3=—r p=—r 


(33) 

(34) 


We concentrate on the first inequality (33). When Xi is 0, the inequality is trivially fulhlled. 


Thus, we can restrict ourselves to the case when Xi > 0, and hence Xi = Xi- If we set e = 0 and 


replace (x, y) by {x,y) everywhere, the difference between the two sides of (33) is the quantity 


Qx in Lemma [T6| and hence it is bounded by K. Going back from {x,y) to (x,y) cannot make 
the right-hand side smaller. Thus we are done if we prove that the “slack term” term e ■ Xi is 
at least K. This is true by construction, since the non-zero values of Xi are at least K/e. The 


other inequality (34) follows similarly. 


This concludes the proof of the lower bound. 


□ 
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The theorem can be extended to more than two coupled recursive sequences. Then we need 
a separate parameter 5 for each function in (28 29) These parameters must be determined from 
a system of equations, and the no-drift condition ensures that this system has a solution. 

The technical condition of Theorems and that certain coefficients are non-negative 
have the purpose to exclude periodicity and can be replaced by weaker conditions. 


6.4 Asymptotic growth constants 

We apply Theorem |13| to the recursion describing the r-chain with corners. It is straightforward 
to compute the 2x2 condensed coefficient matrix (22) with a computer by accumulating all terms 


derived in Section 6.1 and to compute its dominant eigenvalue. Since n = rk 1, the growth 
constant Tj. in terms of n is r-th root of this eigenvector. We observe the same phenomenon 
as for chains without corners, see the right-most column of Table The values increase to 
some maximum, and then the taper off and converge to 3 as r increases further. The first two 
entries in the table reproduce the results for the double-chain (the condensed coefficient matrix 
is ( 22 ) which gives Ti = 3) and the double zigzag chain (the condensed coefficient matrix is 
(ID which gives T 2 ~ 3.0532). We observe that the maximum is achieved for the 8-chain with 
corners (r = 8). 

To establish this bound rigorously as a lower bound, we have to check the conditions of 
Theorem 13 It is easy to check that the coefficients are indeed positive for 


/? = — 1, OTTi The condensed coefficient matrix (22) is 


A^c ^CF\ ^ /2835 2619\ 
AFC aFFJ ~ 1^6022 5504 ) ' 


(35) 


Its dominant eigenvalue is M = (8389 -|- \/69945633 )/2 « 8376.175, with corresponding left 
(unnormalized) eigenvector {px,Py) = (6022, M — 2885) and right eigenvector = 

(2619, M — 2885)"'". The matrix of total group-to-group jump sizes is 

(DCC dCF\ _ /-2619 0 \ 

\dFC dFFJ - 1^-2619 2619y ’ 


Weighting these numbers with the eigenvectors and summing them up (23) yields that the 
weighted total jump size D is zero. The conditions of Theorem 13 are thus fulfilled. 

An intuitive explanation of the equality D = 0 might be as follows. The recursion between 
the two vectors and is not symmetric, as witnessed, for example, by the non-symmetric 
condensed matrix (35). This asymmetry comes from the arbitrary decision to cut the construc¬ 
tion to the right of each corner point. However, on the whole, this irregularity should not cause 
a systematic “drift” in the recursion, which would favor a tendency towards larger or smaller 
numbers i of unhnished runners crossing the cut. Thus, it is not surprising that D = 0. We 
expect that D = 0 should hold for all r, but we have only checked it numerically for small values 
of r, and we have established it rigorously only for the concrete case r = 8. 

the sequence Fq grows at most like M^ and at least like (M — e)^, for any 


By Theorem 


13 


e > 0. Since n = 8A: -|- 1, the growth constant in terms of n is Tg = v^M ps 3.093005695. 


Corollary 17. The 8-chains with corners have 0{Tg) and ^((Tg — e)^) down-free matchings, 
for every e > 0. □ 

This implies Theorem with the help of Theorem 

Numerical data suggest the more precise estimate Fq = M^+ ci/k + 0{l/k‘^)) 
with Co ~ 0.1321 and ci ~ —0.102. This has been computed by Moritz Firsching (personal 
communication) by interpolation from the elements Fq^^, Fq^^, Fq ^'^,... ,Fq^^^, assuming that 
the sequence has the asymptotic form Fq = C^fk°‘{co c/k C 2 /k‘^ This method has 
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predicted more than 300 correct decimal digits of C = M, and it gives a = 3j2 also to a precision 
of more than 300 digits. By comparison, for the sequence ak of down-free matching numbers 
of the zigzag-chain (Section]^, for which the explicit generating function and hence the form 
of the asymptotic growth is known, the same method gave estimates that were accurate also to 
more than 300 digits, both regarding the growth constant l//i = (\/93 -|- 9)/2 and the power 
a = 3/2 of the polynomial factor. 

The asymptotic growth of the form = cqM^/ l?l‘^(\ P o(l)) is not unexpected; it is in 
accordance with the behaviour of r-chains without runners, which has been derived in the proof 
of Proposition]^ (Section |5.4[) by lattice path method [3l Theorem 3]. 


7 Concluding remarks 

7.1 Table of results for pm, dfm and am 

In Table we summarize asymptotic bounds on different structures for three kinds of matchings 
considered in this paper - pm, dfm and am. Some of them do not follow from results proven 
or mentioned in this paper, and we explain them below. First we want to point out some 
observations that can be seen in the table. 

Obviously, pm(X,^) < dfm(X,j) < am(X„), but is dfm(A'„) more likely to behave similarly to 
'pm[Xn) or to am(Xn)? Table shows that different possibilities exist. For a downward chain 
SCn, every matching is down-free and thus dfm(SC,i) = am(Ai,i), but for an upward chain dfm 
is equal, up to a polynomial factor, to the lower bound. For SZZCn, the three growth constants 
are all different, but the intermediate basis for dfm is closer to the upper bound. However for 
r-chains withont corners, as r grows, the growth rate for pm and dfm tends to the same constant, 
3, from below and from above respectively; whereas that for am tends to 4. 



pm(A„) 

dfm(A„) 

am(Xn) 

SCn upside down 

Cn/2 = &*m 

(i„/2i) = e'( 2 ’‘) 

Mn = 0*(3”) 

SC,. 

Cn/2 = e*{2n 

Mn = 0*(3”) 

Mn = 0*(3”) 

SZZCn 

0* (2.1974”) 

0*(3.O532”) 

0*(3.1O22”) 

CH*(ll,n/ll) 

0*(2.5517”) 

0* (3.0840”) 

0*(3.4614”) 

CH*(r,n/r), r —)• 00 

0*(a”), a / 3 

0*(/3”),/3\3 

0*(7”), 7/4 

CH( 8 ,(n- l)/ 8 ) 


0* (3.0930”) 


CH(r, (n — l)/r), r —)> 00 


0*(<5”), (5\3 ? 


DC„ 

0*(3”) 

? 

0*(4”) 


Table 2: pm, dfm, am for several structures. 


Now we describe the entries of the table. The first two lines are classical results, except for 
the formula dfm = ™ upward chain, which has been proved in Proposition 


10 


The estimate pm(SZZCn) = 0*(2.1974”) from [T] was mentioned in Section]^ Actually, it 
was the fact that pm increases from SC to SZZC which initially prompted us to try whether the 
old record of the double structure DC could be beaten by the corresponding double structure 
DZZC. The formula dfm(SZZCn) = 0*(3.0532”) is the main result of Section]^ The estimate 
am(SZZC,i) = 0*(3.1022”) can be derived in a similar way, by adding an appropriate term to the 
recursion i) for Ofc: the only difference is that when Pi is matched to P 3 , the point P 2 can be free. 
The singularity closest to 0 of the resulting generating functions occurs now in (\/l05 — 9)/12, 
one of the roots of 1 — 9x — 6 x^. Thus, in this case the base is 3.1022. 
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For r-chains without corners, CH*(r,/c), the growth rate for dfm has been determined in 
Section ^ and, as was discussed in Section [5^ it converges to 3 from above as r —)• oo. The 
other entries in the line for CH* can be obtained by modifying the analysis of Section |5.3[ we 


only need to replace appropriately in the formula for zi in Proposition 


10 


the factor 


.L(r-i)/2j 




representing the number of down-free matching on an arc of r — i points. For pm, we have to 
replace it by the Catalan number C(^j.-i )/2 when r — i is even and by 0 when r — z is odd; for 
am, we replace it by the Motzkin number Mr-i- The row sums of the recursion matrix can be 
obtained by plugging these modified expressions for z} into (12). For pm, the resulting sequence 
of row sums is the sequence A189912 from [8], and for am, it is the sequence A077587 (we omit 
the proofs). From the asymptotic behavior of these sequences it follows that their r-th roots, 
which are the growth constants, converge to 3 and 4 from below. 

The growth rate for dfm for r-chains with corners, CH(r, A:), was treated in Section]^ Em¬ 
pirically, they seem to be better than r-chains without corners. The monotone convergence to 3 
from above is not proved. It seems plausible that the difference between r-chains with corners 
and r-chains without corners should become negligible as r —)> oo, and therefore the growth con¬ 
stant should converge to the same constant 3. That the convergence should be monotonically 
decreasing is only based on the empirical observation from Table We have not extended the 
analysis to pm and am, although this would be feasible with some effort. We expect that the 
results would be the same as for r-chains without corners. 

The formula pm(DCn) = 0*(3"^) is the classical result of Garcia, Noy, and Tejel [6], in 
accordance with dfm(SCn) = 0*(3"^) from the first line. The estimate am(DCn) = 0*(4”) is due 
to Sharir and Welzl [9], and it is currently the best lower bound on the maximum number of 
am. The growth of dfm(DCn) remains unknown, but it is hi*(3”) and 0*(4’^). 


We see no reason to think that our best construction CH(8, h) is optimal in the sense that it 
has the maximal possible dfm and/or that the corresponding double construction has the max¬ 
imal possible pm. Sets with asymptotically higher bounds may very well be more complicated 
— both in terms of their description and their analysis. An obvious continuation from single 
chains to r-chains would be to insert a third level of downward arcs between the vertices of 
r-chains, possibly continuing towards a fractal-like pattern. We have not attempted to analyze 
these structures. 


7.2 Summary and Outlook 

We have found new constructions of point sets with a larger number of perfect matchings than 
previously known. More importantly, we show that, like for triangulations, the true bound 
for perfect matchings is not given by the double chain. For the analysis of these sets, the 
notion of down-free matchings was crucial. It allowed us to concentrate on one half of a double¬ 
construction. 

We have also shown that methods from analytic combinatorics are useful for counting prob¬ 
lems for geometric plane graphs. However, the results from analytic combinatorics that we are 
aware of cannot be readily applied for r-chains with corners. In this case, the analysis leads to 
coupled recursions involving two sets of variables. For these recursions, we had to develop our 
own methods. These somewhat pedestrian methods give the growth rate only up to an arbi¬ 
trarily small error e. We hope that that the methods of analytic combinatorics will be further 
developed to encompass such cases as well. 
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